On the Decoherence of Primordial Fluctuations During Inflation 
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We study the environment-induced decoherence of cosmological perturbations in an inflationary 
background. Sphtting our spectrum of perturbations into two distinct sets characterized by their 
wavelengths (super and sub-Ifubble), we identify the long wavelength modes with our system and 
the remainder with an environment. We examine the effects of the interactions between our system 
and the environment. This interaction causes the long-wavelength modes to decohere for realistic 
values of the coupling and we conclude that interactions due to backreaction are more than sufficient 
to decohere the system within 60 e-foldings of inflation. This is shown explicitly by obtaining an 
analytic solution to a master equation detailing the evolution of the density matrix of the system. 
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I. INTRODUCTION 

As is well known, temperature fluctuations in the CMB 
and the inhomogeneities that seed structure formation 
in the universe share a common origin. Both are a re- 
sult of the scalar metric perturbations produced during 
inflation. However, these perturbations are of a purely 
quantum mechanical nature while no cosmological sys- 
tems of interest (CMB anisotropics, clusters etc.) display 
any quantal signatures. Presumably, for this to be the 
case, the primordial density perturbations underwent a 
quantum-to-classical transition some time between gen- 
eration during inflation and recombination, when struc- 
ture flrst became apparent. 

Decoherence is a much studied process (see Q for a 
comprehensive review). Although not all conceptual is- 
sues have been resolved, it is understood that it can oc- 
cur whenever a quantum system interacts with an "en- 
vironment" . In other words, this effect can be said to 
pervade open systems due to the difficulty of creating a 
truly closed, macroscopic quantum system. Along with 
its ubiquity, it is also known to be a practically irre- 
versible process, since the loss of quantum correlations 
in the system is accompanied by an increase in entropy. 

Early studies of the classicalization of primordial per- 
turbations focussed on intrinsic properties of the system 
(see, for example 2], 3]). This was made possible by the 
application of ideas of quantum optics to the theory of 
cosmological perturbations. Primordial density fluctua- 
tions (the scalars as well as the tensors) evolve into a 
peculiar quantum state - a squeezed vacuum state 
By studying the large squeezing limit of these states, it 
was found that quantum perturbations become indistin- 
guishable from a classical stochastic process. In other 
words, quantum expectation values in a highly squeezed 
state are identical to classical averages calculated from 
a stochastic distribution, up to corrections which vanish 
in the limit of infinite squeezing. The authors of Q refer 
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to this as "decoherence without decoherence" while [6j 
endows the phenomenon with the more technical epithet 
" quantum non-demolition measurement" . We emphasize 
that these works focussed on the classical properties of 
the states and not on the coherence properties of the sys- 
tem. 

As is well understood, in order to study true classi- 
calization, one must consider two distinct aspects of a 
system. First the quantum states must evolve, in some 
limit, into a set of states analogous to classical configura- 
tions. The second is that these resultant states interfere 
with each other in a negligible fashion. This last property 
constitutes decoherence and is equivalent to the vanish- 
ing of the off-diagonal elements of the density matrix. 

A truly closed gravitational system is a practical im- 
possibility (unless one considers the totality of the uni- 
verse to constitute the system as in, for example, quan- 
tum cosmology). Since the gravitational interaction has 
infinite range and couples to all sources of energy, inter- 
actions with some sort of environment are an inevitabil- 
ity. As such, environmentally induced decoherence must 
also be present and would play an important role in the 
classicalization of primordial density fluctuations. 

The purpose of the present article is to determine pre- 
cisely the effects by the "inflationary environment" (we 
will elucidate this notion below) on cosmological pertur- 
bations and to study the resultant decoherence. Other 
authors have also examined this problem (see, for exam- 
ple 0],[3,[1|,|3jE3) " however, we are the first to present 
an exact analytic expression for the density matrix with 
a realistic environment-system interaction. 

The paper is organized as follows: in the next section, 
we review some basic properties of decoherence of which 
we will make use. After reviewing the quantum theory of 
cosmological perturbations in section III, we make clear 
our concept of the environment and motivate some real- 
istic interactions in section IV. Subsequently, we develop 
necessary formalism which, in section VI, we make use 
of to demonstrate the classical nature of the system and 
calculate the decoherence time scale. 
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II. DECOHERENCE 

In the present section, we intend to present an ex- 
tremely (but, hopefully, not exceedingly) terse account 
of the theory of decoherence. The physics of classicaliza- 
tion is elegant and subtle and a thorough exposition of its 
finer points would bring us too far afield from the purpose 
of this article. We confine our attention solely to the car- 
dinal features and disregard any peripheral aspects. The 
reader unsatisfied by our presentation is encouraged to 
consult any of a number of excellent reviews of which we 
mention but a few [lll| . [l^ . [H . 

From an operational perspective, the process of de- 
coherence usually refers to the disappearance of off- 
diagonal elements of the density matrix. These elements 
(phase relations) represent the interference of states in- 
herent in any quantum system. Evidently, their disap- 
pearance is an integral part of a quantum-to-classical 
transition. 

Having mathematically defined decoherence, we now 
turn to the physical processes responsible for it. At the 
heart lies the concept of the open system and the near 
impossibility of forming a macroscopic closed state. Vir- 
tually all realistic systems must interact with an envi- 
ronment of some sort where, by environment, we refer to 
degrees of freedom which interact with degrees of free- 
dom in our system but which are not witnessed by some 
observer intent only on the evolution of the system. This 
leads to the first important characteristic of decoherence 
- its ubiquity. 

Next, we come upon the concept of entangled states. 
Initially, if we disregard all correlations between system 
and environment, our composite wave function (system 
-|- environment) can be expressed as the outer product 
of the system and environment states (more generally, 
it will be the outer product of ensembles of states, as 
is the case when one makes use of density matrices). 
Though initially factorizable, interactions between the 
environment-system pair rapidly change this: the total 
state evolves from the form 

i j 

to 

i,j 

where ^ represents an entangled state and, as such, is 
non-factorizable in this basis. Entanglement is key to 
the whole process for the following reason - an entangled 
state produces a density matrix which is non-factorizable. 
The operational equivalent of an observer ignoring the 
environmental degrees of freedom is to trace out (partial 
trace) these degrees of freedom. Due to the orthogonal- 
ity of the environment states, the observer is left with a 
density matrix which diagonalizes as the states entangle 



(the fact that the decoherence rate is related to the rate 
of entanglement has been used to estimate decoherence 
times. See, for example, @,0|). An interesting prop- 
erty of classicalization follows from this - the interference 
terms are still present, but are unobservable by a "lo- 
cal" observer (local in the sense that he only observes 
the system). 

These "hidden" interference terms lead us to our next 
point. By tracing out the environmental degrees of free- 
dom, an observer throws away all the correlation terms, 
leading to a decrease in the amount of information avail- 
able in the system - hence, this leads to an increase in the 
entropy from which we can conclude that decoherence is 
a practically irreversible process. 

The system being decohered, it can only be found in 
a much smaller subset of the states that were previously 
allowed - this is what prevents us, in part, from seeing 
"Schroedinger's Cat" states at a macroscopic level. The 
states that diagonalize the density matrix of the system 
are referred to as pointer states and these states 

remain in the subset of physical states after decoherence. 
If the evolution of the system is dominated by the self- 
Hamiltonian of the system, the pointer basis is composed 
of the eigenstates of the self-Hamiltonian while, if the 
interaction dominates, the eigenstates of the interaction 
form the basis Pointer states are also those states 

for which the production of entro py during decoherence 
is minimized (predictability sieve) 16]. 

Finally, we conclude with a heuristic view of decoher- 
ence. Neglecting certain interacting degrees of freedom in 
a theory will generally lead to an apparent loss of unitar- 
ity. Thus, one should expect a flow of probability out of 
the system which, in turn, manifests itself as a vanishing 
of certain elements of the density matrix. 

III. QUANTUM PERTURBATIONS IN AN 
INFLATIONARY UNIVERSE 

A. The Action for Quantum Perturbations 

We provide in this section an overview of the quantum 
theory of cosmological perturbations in an inflationary 
background. For a more in-depth treatment, the reader 
is referred to or ITij. 

The classical action for an inflationary model is given 
by (in this and in what follows, we set G = fi. = 1) 

S ^ J d^x^i^R+^d^cl^d'^cl.-Vm. (3) 

If the potential V{(j)) for the matter scalar field (j) is 
sufficiently flat and if, in addition, initial conditions are 
chosen for which the kinetic and spatial gradient terms 
in the energy density are negligible, this action leads to 
a period of inflation during which the space-time back- 
ground is close to de Sitter 

ds^ = {-f{-drj^ + {dxT), (4) 
V 
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where rj is conformal time. 

During the course of inflation, any pre-existing classical 
fluctuations are diluted exponentially. However, quan- 
tum fluctuations are present at all times in the vacuum 
state of the matter and metric fluctuations about the 
classical background space-time. Their wavelengths are 
stretched exponentially, become larger than the Hubble 
radius H~^{t) and re-enter the Hubble radius after infla- 
tion ends. These fluctuations are hypothesized to be the 
source of the currently observed density inhomogeneities 
and microwave background anisotropics. In order for this 
hypothesis to be correct, the fluctuations must decohere. 

The quantum theory of linear fluctuations about a clas- 
sical background space-time is a well-established subject 
(see e.g. the reviews ^3 If the matter has no 

anisotropic stress (which is the case if matter is described 
by a collection of scalar flelds) , then a gauge (coordinate 
system) can be chosen in which the metric including its 
(scalar metric) fluctuations {^) can be written as 

ds^ = (-)2(-(l + 2V'(a:,77))V + (l-2^(x,7?))(dx^)2), 

. . . . . 

and the matter including its perturbation [6(j)) is 

^^^ + S(t>{x,7i). (6) 

The quantum theory of cosmological perturbations is 
based on the canonical quantization of the metric and 
matter fluctuations about the classical background given 
by a(ri) and Since the metric and matter fluctua- 

tions are coupled via the Einstein constraint equations, 
the scalar metric fluctuations contain only one indepen- 
dent degree of freedom. To identify this degree of free- 
dom, we expand the action (0) to second order in Scj) and 
■0, and combine the terms by making use of the so-called 
Mukhanov variable [HE^I 

V = a{ri) [5(j)+ ^V], (7) 

in terms of which the perturbed action S2 takes on a 
canonical form (the kinetic term is canonical) and the 
perturbations can hence readily be quantized: 

S2^\jd^x[v'^-v,,v^, + ^v\ (8) 

where z = and a prime indicates a derivative with 
respect to 77. This action contains no interaction terms: 
it represents the evolution of a free scalar fleld with a 
time-dependent square mass 

7" 

m2=-— , (9) 

propagating in a flat, static spacetime. This action leads 
directly to a well-deflned quantum theory via the canon- 
ical commutation relations. 



The Hamiltonian corresponding to the above action ^2 
can be written down in second quantized form: 

H = j £k[k{a\a^ + aKa_^ + l)-i^{aj: a_^-h.c.)]. 

(10) 

The first term in the brackets represents back-to-back 
harmonic oscillators, in phase such that the system has 
no net momentum. The second term leads to the "squeez- 
ing" of the oscillators on scales larger than the Hubble 
radius H~^{t) (on these scales the second term in (|10f) 
dominates over the spatial gradient terms coming from 
the first term in the equation of motion for v). On these 
scales, the squeezing leads to an increase in the mode 
amplitude 

^kijl) ~ z{ii) - a{ri) , (11) 

where the second proportionality holds if the equation 
of state of the background geometry does not change in 
time. We take this to be the case in our subsequent 
analysis. 

B. Properties of Squeezed States 

There exists an extensive literature on squeezed states. 
We refer the reader to [23| and |^] for the mathematical 
properties of squeezed states. For their physical interest, 
we direct the reader to [23l |. 

The evolution of a state of a system governed by the 
Hamiltonian H1Q|I can be described by the following evo- 
lution operator: 

U = S{ru,Vk)R{ek), (12) 

where 

5,.(^) = exp[^(e-2^'^^^('')a_,.a,~ - h.c.)l (13) 

and 

i?(6lfc) = exp[-i6lfe(4afc + aLfca_fc)], (14) 

where S{rk,(pk) is the two-mode squeeze operator, R{Ok) 
is the rotation operator, and the real number r is known 
as the squeeze factor. The rotation operator and the 
phase 9k play no important role in what follows hence 
we ignore them from now on. 

The action of the squeezing operator on the vacuum 
results in squeezed vacuum states 

Skim ^ \k) = (15) 
cosh(^.(,)) (-^'"^'°^''^ tanh(r,(r7)))" |n,fc;n,-/c > . 

The behaviour of the squeezing parameter is com- 
pletely determined by the background geometry. The 
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evolution of the squeezing parameters is typically very 
complicated, but an exact solution is known in the case 
of a de Sitter background : 



= -J-iarctan(^), 



9fe = A:ry + tan-i( — ), 



(16) 



(17) 



(18) 



where the vacuum state being operated upon corresponds 
to the Bunch-Davies vacuum. 

The squeezing operator has the property of being uni- 
tary so that 



1. 



(19) 



Although squeezed states do not provide a basis (as they 
are overcomplete) , they do form an orthogonal set of 
states: 



(20) 



This follows from the properties of many particle states. 

An important property of squeezed states of which we 
will make use is the fact that the number of particles in 
such a state can be expressed entirely in terms of the 
squeezing parameter via 




FIG. 1: The Penrose diagram for de Sitter space 
in planar coordinates. Note that these coordi- 
nates only cover half the spacetime. Blue lines 
indicate lines of constant t, red lines constant r, 
and the solid black line represents the horizon. 



{k\Nk\k) = sinh2(rfc), 



(21) 



where iV/jis the number operator for the k-mode. Phys- 
ically, squeezed states represent states which have min- 
imal uncertainty in one variable (high squeezing) of a 
pair of canonically conjugate variables - the uncertainty 
in the other is fixed by the requirement that the state 
saturates the Heisenberg uncertainty bound. For those 
states of cosmological interest, we take the squeezing to 
be in momentum. 

For our application, the squeezing parameter will be 
quite large. As shown in p^ . 



this section, we show that de Sitter space naturally pro- 
vides us with a hidden sector and that the borderline 
between the visible and invisible in our theory is natu- 
rally given by the Hubble scale. 

Although de Sitter space is geodesically complete, a 
geodesic observer will be subject to the effects from both 
a particle horizon and an event horizon H^,!!^. That 
the latter constitutes a true event horizon can best be 
seen by examining the behaviour of null geodesies in 
Painleve-de Sitter coordinates (see, for example [l^l), 
which remain finite across the horizon, in contrast to 
static coordinates. Specifically, we have 



a[ti) 



(22) 



where a{ti) (0(^2)) is the scale factor at first (second) 
Hubble crossing. For current cosmological scales, k, 
102. 



C. The Hidden Sector 

An essential ingredient in the theory of decoherence 
is the presence of unobserved, "hidden" degrees of free- 
dom: their interaction with our system degrees of free- 
dom causes the delocalization of the phase relations. In 



2jdtdr 



'dn^. (23) 



Here, 1(= 1/H) denotes the de Sitter radius. Clearly, 
setting r — I causes our timelike coordinate to become 
spacelike (the characteristic feature of an event horizon) . 
Timelike observers that cross from r — \e\ to r + \e\ find 
themselves incapable of getting back, trapped outside of 
a sphere of radius 1. 

Now, if one transforms to the coordinates typically 
used when discussing inflation (the so-called planar 
coordinates) and examines the behaviour of timelike 
geodesies, one finds that all timelike observers originating 
within the horizon must eventually cross. 
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The zero-point fluctuations induced by the horizon Ijg 
can be thought of as the seeds for metric perturbations 
, [s^l • HeuristicaUy, the horizon can be thought of as 
a source of thermal radiation with a temperature H/2'it 
(in complete analogy with the black hole case) . This ra- 
diation then produces gravitational metric perturbations, 
with the same spectrum, which are stretched out by sub- 
sequent cosmological evolution and ultimately lead to the 
formation of structure in the post-inflationary universe. 



Note, however, that this naive picture is not quite cor- 
rect - the equation of state of the produced radiation 
is not thermal [sj, and including the effects of gravita- 
tional back-reaction leads to corrections to the thermal 
spectrum (this is also true in the black hole case [3^). 
However, our ensuing discussion in no way relies on strict 
thermality. 



We consider our observer to be to the left of the hori- 
zon in fig.l. In accord with our discussion above, we 
take our radiation to be produced at the horizon with 
a continuous distribution such that a non- vanishing sub- 
set of our modes have wavelengths less than 1 (or H~^). 
It follows that our observer in planar coordinates, due 
to the event horizon, will be prevented from observing 
certain radiation modes. We conclude that those modes 
which are unobservable are those associated with physical 
wavelengths less than the horizon scale. Of course, gravi- 
tational redshifting will cause these modes to stretch and 
eventually cross the horizon. The point is that particle 
production is a continuous process and we expect that, 
at all times, a certain set of modes will be unobservable, 
and these modes will be associated with physical wave- 
lengths less than . As a result of this, decoherence 
is an inevitability and we define our environment to be 
a set of modes whose physical momenta are greater than 
the Hubble scale. 



Having identified the modes of the theory which we 
must trace out, we ask what happens if we trace out 
additional modes. For example, if an observer was only 
interested in very low energy modes (fc <C H) he could 
ignore (or trace out) modes with {k < H, but not k <C 
H) - surely this would provide an additional source of 
decoherence as it increases the environment. However, 
compare this to the case of an observer who is interested 
in all super-Hubble modes. The second observer would 
see less decoherence than the first. Decoherence is, after 
all, an observer dependent effect - an observer who could 
monitor every degree of freedom in the universe wouldn't 
expect to see any decoherence. However, our goal is to 
determine a lower bound on the amount of decoherence 
as measured by any observer in the " out" region of our 
Penrose diagram. In this case, we trace out only those 
modes which we must (i.e. all modes on sub-horizon 
scales) and take our system to be composed of the rest. 



IV. INTERACTIONS WITH THE 
ENVIRONMENT 

Key to our investigation of decoherence is the notion 
of the environment. Such an environment can take on 
many different guises. As was stated above, we define 
ours in the following fashion: expanding the background 
fields (gravity and the infiaton) in terms of fiuctuations, 
we identify our environment with the fluctuations whose 
wavelengths are less than some cut-off, while our system 
consists of those wavelengths greater than this cutoff. As 
explained above, since we are operating in a de Sitter 
background, the natural scale to pick for the cutoff is the 
Hubble scale. 

In order to determine the precise form of interactions 
inherent to a system of cosmological perturbations, we 
expand ((SJl to the next order (recall that expanding to 
second order is what led to a free field theory) in the 
fluctuations, and express the result in terms of v{x,r]). 
Interactions can either be purely gravitational in nature 
(backreaction), or they can arise in the matter sector 
through V{(j)), the inflaton potential. 



A. Gravitational Backreaction 

To focus on the interactions due to gravitational back- 
reaction, we must expand the gravitational action to 
third order in the amplitude of the perturbations and 
write down the potential in terms of the Mukhanov vari- 
able V. Expanding to higher order simply introduces 
more complicated interactions. For our purposes, we re- 
strict our attention to the simplest terms that arise. 

In the case where the metric, including its fluctuation 
field "0, is given by 



a2(7^)[-(l + 2'ilj)dT]^ + {1- 2'4j){dxy 



(24) 



we can expand the Ricci scalar in powers of tp to obtain 



R = 



6air,y 



6«(1 



2tP + 4tP'^ -8^p^...) 

(25) 

(where terms with derivatives either temporal or spa- 
tial of the ijj have been ignored as they are sub-dominant) 
from which we can extract our term of interest, R^^\ 



(26) 



which is the leading order gravitational self-interaction 
term. Recalling the definition ^ of the Mukhanov vari- 
able in a slow-roll infiationary background, our potential, 
expressed in terms of v, becomes (neglecting S(j) when 
substituting v for ^ and we use the fact that, for our 
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inflationary background, a(?7) = l/{Hri)) 



V = 



IGttM 
1 



Ml J 



\/27r 



3 a(77) 3 



so that 



V 



'Mpi (2e)3/2 



TT a'^(?7) {(j))'a 



with 



A 



3 1 



Mpi (2e)3/2 

and where we've used the slow roll conditions 



and 



2 ' 



(27) 



(28) 



(29) 



-a{ri) = a(77)Ao, (30) 



ff2 = V{(j))/{-iMli), 3H4) = -t/', (31) 



(32) 



is one of the slow-roll parameters. Our dimensionful cou- 
pling is explicitly time-dependent - this is to be expected 
since it is associated with a fixed physical scale and our 
theory ((SJ is written entirely in terms of co-moving quan- 
tities. 



Mukhanov variable (and this time neglecting -0 in the 
process of substitution) , we have 



I d?xAa'^{'n)H(t>{-f = I d''xa{r])Ati(l3v'', (35) 



so that 



A — A^(j)a{rii). 



(36) 



How do the coupling strengths of the two potentials 
compare? Taking the ratio of the two, we find 



A 



inf 



4^0a(?7) 



37r 1 
V2Mpi (2e)3/2 



(37) 
(38) 



Since the observationally allowed value for (p at times 
when fluctuations relevant to current observations are 
generated is of the order 10~^Afp/, we find that the grav- 
itational coupling could conceivably dominate depending 
on the value of e. Since we are only interested in ob- 
taining a lower bound on the decoherence rate, and due 
to the fact that the exact form of the inflaton potential 
(along with the initial conditions that determine e) is 
model dependent, we consider gravitational backreaction 
to be the main source of decoherence in what follows. 
Nonetheless, the above demonstrates that inflaton inter- 
actions have the potential to be important. 

We couple our system to the environment by writing 



B. Inflaton Interactions 

In addition to the gravitational backreaction terms, 
there are also interactions due to non-linearities in the 
matter evolution equation. Consider a model of chaotic 
inflation with a potential of the form 




where /i is a dimensionless coupling constant. The per- 
turbations produced during inflation are joint matter and 
metric fluctuations. The matter part of the fluctuation, 
denoted by 5(j), give rise to a cubic term in the interaction 
potential of the form 

y - y d?x4,yf^^i(l){5(j)f, (34) 

where, in the case of slow-roll inflation, we can treat (j) 
as a constant. Now, writing the potential in terms of the 



V = j d^xXv^ = j d^xXv^ip, (39) 

where v now refers only to the expansion of the 
Mukhanov variable in momenta greater than some cut- 
off and if is the same field but expanded in terms of the 
environment modes. 



V. THE DENSITY MATRIX 

Having determined a candidate interaction between 
our system and the environment, we now face the task of 
deriving an appropriate master equation in order to de- 
termine the time dependence of our density matrix. Sev- 
eral approaches exist (for example, H^jfl^) which have 
been used by a number of authors - rather, we follow the 
method of (35j which we now review. 

We assume that our system of interest is weakly in- 
teracting with some environment. The Von Neumann 
equation for the full density matrix (p) reads (note that 
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we make use of conformal time. This is due to the fact 
that our action (jHJ is expressed in terms of conformal 
time) 



^ = -^[H.p\, (40) 

where H is the total Hamiltonian of the system and can 
be written as 



We find that (see the appendix) 

TrB{V{xi,'qi)V{x2,Tl2)p) = ^'^ ^^^^ Km ^ m)5{xi - X2). 

(46) 

In light of the fact that this equation was derived in 
the limit of small time intervals, we can approximate the 
integral in (|45|l by the product of the integrand with the 
time interval, rj. Bringing (po) to the left-hand side and 
dividing both sides by time allows us to write 



H = Ho + V, 



(41) 



where Hq is the self-Hamiltonian and V couples the sys- 
tem to the environment. Note that (p) denotes the full 
density matrix for the system and the environment. 

Switching to the interaction representation H40|l takes 
on the form 



where 



TTJ -1 

— = -i[V,p], 

arj 



p = exp{iHor]) p e:cp{-iHor]), 



(42) 



(43) 



with a similar expression for V. 

A perturbative solution of is found to be given by 
the following: 



Po-i dr [V{t),po] 



(44) 



Jo Jo 



Our ultimate goal is to derive an equation of motion for 
the reduced density matrix {pA = Trsp, where A denotes 
the system quantities while B refers to the environment. 
We use this notation throughout the rest of the paper) . 
To this end, we trace out the environmental degree of 
freedoms to obtain 



PA 



P^ 



dT2 / dnTrB[V{T2),[ViTi),po]] 
Jo 



(45) 

Note that the first order term has vanished - this is due 
to the specific form of our system-environment coupling. 
Had we used a potential in which an even power of the 
environment field had appeared, we would have obtained 
a non- vanishing contribution at this order. Had this been 
the case, the first order term could have been neglected 
on the grounds that it would lead to unitary evolution of 
the system - since our goal is to study the decoherence of 
the system (a non- unitary process), we can safely ignore 
such terms. 



p- Po _ dp_ 
rj dr] ' 



(47) 



in the limit of small rj. 

As the initial time (77 — 0) is arbitrary, we conclude 
that our equation for the reduced density matrix may be 
written as (in terms of physical time) 



dm 

dt 



-a{t) 



d^x[v^,y,p{t)]h (48) 



where the details have been relegated to the appendix. 
V is a normalization volume, Ti. — a(t)H, with H the 
physical Hubble scale, and we note that in order to obtain 
the condition H46|l , it was necessary to eliminate non-local 
terms by coarse-graining over scales of order Ti. in both 
time and space. 

The differential equation 1)48(1 is the master equation 
for our system. In order to proceed, we obtain a matrix 
representation in the basis of squeezed states. Again, we 
point out that these do not form a true basis for the 
Hilbert space (note, however that the use of an overcom- 
plete basis poses no difficulties when it comes to obtain- 
ing representations of the density matrix |llj|'). However, 
in the limit of large squeezing, squeezed states become 
orthogonal to other states in the system. Since squeezed 
states are the "natural" states of the system, we view all 
other states as being spurious and truncate our Hilbert 
space so that it contains only the former. Furthermore, 
as our interactions are small compared to ()10ll . we iden- 
tify the squeezed states as our pointer basis 36] . 

Finding a matrix representation of ea. (|48|l is a rela- 
tively simple affair - due to the nature of the squeezed 
states, the expectation value operator w^",n G Z must 
be diagonal in this (discrete) basis of states. This, along 
with the identities 



S{rk,ipk)a±kS^{rk,ipk) = a±fc cosh(rfc) 



sinh(rfc). 



and 



5* ^{rk,ipk) = S''{rk,ipk) = S{~rk,ipk), 



(49) 



(50) 



renders the calculation relatively straightforward. Note 
that {k\Nk\k) — sinh^(rfe), where Nk is the number op- 
erator 21]. With this in mind, we find that (|48|l reduces 
to 
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dt 



, ,128^2 ^2(^) sinh*(r,) 
sinh^ (rj ) sinh^ {rj ) 



sinh* {r.j ) 



(51) 



where, for simplicity, we've replace the cosh^(r) terms 
with sinh^ (r) since we are interested in the limit of large 
r. ^ 

The combination ^"'M''') ^ ni{t) = a^{t)ni{0) is to 
be interpreted as the particle density, a quantity which 
is finite in the thermodynamic limit. Clearly, the deco- 
herence rate increases as the difference between the two 
momenta increases. For this reason, we take our states of 
interest to have approximately the same momenta, and 
the above reduces to 



dt ^ ' 



A2nf(0) (fc.~%)^ 



i?5 



Pv, (52) 



in terms of physical time and co-moving momenta and 
volume. 

A few things are immediately obvious: 

1) The diagonal elements suffer no loss of coherence. 
This actually could have been surmised much earlier from 
eq. (|45|l by noticing that the trace over the system degrees 
of freedom must vanish. 

2) The rate of decoherence grows extremely rapidly. In 
fact, in order to decohere the system within 60 e- foldings 
(approximately the minimal time permissible for the du- 
ration of inflation), the initial particle density (no) can 
be as low as 10~^^ particles per Hubble volume. 39] 

3) The particular time t = for a pair of modes should 
be taken to correspond to the the time that the shortest 
of the pair (the higher energy mode) crosses the horizon. 

So far, we've argued that a certain sector of the the- 
ory is unobservable (thus justifying a minimal amount 
of tracing), determined an interaction between our vis- 
ible and invisible sectors, and obtained a lower bound 
on the parameters of the theory such that decoherence 
takes place within 60 e-foldings of inflation. The ques- 
tion remains: in a realistic cosmological model, are the 
parameters of the theory such that decoherence can take 
place during the inflationary period, and be caused by the 
leading order gravitational back-reaction term? In other 
words, is the bound we found satisfied in conventional 
models? 

In order to answer that question, we must obtain the 
number density of particles in a typical super-Hubble 
mode at first Hubble crossing. 

Consider the square of the substitution we used to ob- 
tain our potential in terms of the Mukhanov variable: 



To determine the number of particles of the v field in 
terms of physically meaningful quantities, we must first 
quantize the Mukhanov field. However, once the the- 
ory is quantized, the expression H53|l is meaningless - the 
left-hand side is an operator, while the right is a classi- 
cal field. In light of this, we follow the usual route [s^ 
in semi-classical gravity and replace v with it's vacuum 
expectation value: 



iv') = «^(ry)(^)V. 
In the limit of large squeezing, we have that 

1 r d^k 



27r3 / k 



-Nkit), 



(54) 



(55) 



where Nk (t) is the number of particles in the k-mode at 
time t, which scales in time as 



Nk{t) oc a\t), 



(56) 



where we now consider only physical (as opposed to co- 
moving as in the previous discussions) quantities. The 
extra factors of a{t) in the particle number appear be- 
cause we are now considering the red-shifting of the mo- 
menta (see (|16f) '). We expect the spectrum to be ex- 
ponentially suppressed at high (sub-Hubble) momenta: 
therefore, to a good approximation, the integral in H55|) 
can be taken to be over the infrared sector only. Fur- 
thermore, rather than performing the integral over the 
modes, we reparameterize and integrate over the times 
which these particular modes first crossed the horizon. 
In other words, we let 



H 



and 



(57) 



(58) 



where, as above, t — denotes first Hubble crossing for 
a particular mode. We now have. 



H 

dta^{t) = —NH{0)a^itr), 



(59) 

with tj, denoting the time of reheating and where we've 
ignored the time-dependence of the Hubble scale. 

During reheating, the inflaton will undergo periods 
when it's total energy is dominated by it's kinetic term. 
So, during reheating, we can make the substitution (jp ~ 
Pr to obtain 



= a^(^)(^)V. 



(53) 



i?2 



(60) 
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We identify N^H^ — njf{0) with the number of par- 
ticles of momentum H per Hubble voume and taking the 
reheating temperature as H so that pr — i?^. We can now 
make use of the fact that, observationally, ip'^ « 10~^, 
to deduce that n/f (0) « 10^^ particles/Hubble volume. 
This is well above the lower bound we found. In this case, 
we find that the modes will decohere approximately 20 
e-foldings after crossing the horizon. 



VI. CONCLUSION 

In this paper, we have studied the decoherence of cos- 
mological fluctuations during a period of cosmological in- 
flation, taking the effects of squeezing into account. We 
have determined realistic interactions for our system of 
perturbations and have found that, at the same order, 
gravitational interactions and matter (inflaton) interac- 
tions are comparable, depending on the scale of inflation 
and the slow-roll parameter e. Furthermore, we have jus- 
tified the use of Hubble scale as a cutoff. 

Having considered the leading order gravitational cor- 
rection to the action of quantized cosmological pertur- 
bations, we find that super-Hubble modes decohere long 
before the end of inflation. Of course, we have only ob- 
tained a lower bound on the decoherence rate - inter- 
actions more complicated than the ones considered here 
will generally lead to much faster decoherence times |13| . 



VII. APPENDIX: TRACING OUT THE 
ENVIRONMENT 

In this appendix, we explicitly calculate the partial 
trace of ea. (|45|l . 

The expansion of the Mukhanov variable in a spatially 
flat background takes the form 



= (afce-*'="+4e*^-^), (61) 



(27r)3/2 



and we restrict our attention to modes within a sphere 
of radius H in momentum space. Since our calculation 
will be performed in terms of comoving quantities and we 
take our cutoff to correspond to a fixed physical scale, our 
cutoff acquires a time dependence of the form 



H = a{r])H. (62) 
Our normalization conventions are as follows: 



|fc) = v/2^4l0), {k\k') - {2TT)^2EkS'-''\k-k'), 



The identity operator has the form 



d^k 1 



(27r)3 2Ek' 



(64) 



For simplicity, we ignore the effects of squeezing until 
the very last. As our initial conditions, we do not take the 
environment to be in the vacuum - this would be contrary 
to the basic idea of the generation of inhomogeneities. We 
take our states to be 2 particle zero-momentum states. 
Were we to explicitly include the effects of squeezing, we 
would find that our scattering amplitude (i [<?!>"' |j) would 
scales as a" (77). Our approach is as follows: we calculate 
the scattering amplitude for a fixed particle number (2) 
and, at the last step, include the additional factors of 
0(77) in order to embody the effects of particle produc- 
tion (squeezing). Note that we must take into account 
squeezing since ((TB|) tells us that all modes in de Sitter 
space get squeezed. 

We take these states to be populated according to a 
distribution which falls off exponentially in the UV, with 
temperature parameter T = (3^^ — In other words 



Penv = Cexp(-/3i7), 



(65) 



where this H refers to the Hamiltonian. The precise form 
of the distribution is immaterial - after tracing, the only 
information that the systems retains about the environ- 
ment is it's "size" (the cutoff scale). As another simplifi- 
cation, we take the energy of the state to be dominated 
by it's momentum. Due to the nature of squeezing and 
in view of our comments about the distribution, this is 
a perfectly justifiable assumption. C is a normalization 
constant which we determine by the condition that the 
trace of the left hand side of the equation be pays i-e. 



T^^envP 



-^^eni) P 



d^k 1 (fc,-fc|p|fc,-fc) 

{2n)^2Ek {k,-k\k,-k) 



(66) 



C 

^Psys 



dk^exp{-2l3Ek) 

1 



= ^Psvs{^ree-^\l + A^))^Paya. 

Therefore, we set C w leTr^e'^'^/H^. 

The terms on the right hand side will all have the basic 
form (aside from the trace of po, which is the same as the 
above): 



RHS 



d^k 1 {k,-k\v{x)v{x')p\k,-k) 



(2^)3 2Ek 
Stt 



(fc,-fc|fc,-fc) 



(67) 



= (27r)3 5(3)(fc-fc'). 



(63) 



Sim^Jn k{x-x') ^ ' 

where the delta function arises from the normalization 
of the states. Since we are only interested in physics on 
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scales much greater than 7i, we coarse- grain over time 
and use the relation 



Thus, we find that 



(68) 



Note that we identify S{Q) with the volume of space and 
we've included the additional factors of 0(77) as dicussed 
above. 



RHS 



s{o)n' 



■S{v-v') 



H 



k{x — x') 



Again, as our interest lies in scales such that 7i(a; 
xi) ^ 1, we perform the substitution 



^sin[fc(x — a;')] 
k{x — x') 

Finally, we obtain 



) =TT5{Ti.{x-x')). 



RHS 



8772 6{-q-i)5{x-x') 2 



m 



(70) 



a^iv). (71) 
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